Density-functional Theory for / electron Systems: the ct-j Phase Transition in Cerium 
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The isostructural a-y phase transition in cerium is analyzed using density-functional theory with 
different exchange-correlation functionals, in particular the PBEO hybrid functional and the exact- 
exchange plus correlation in the random-phase approximation [(EX+cRPA)@PBE0] approach. We 
show that the Hartree-Fock exchange part of the hybrid functional actuates two distinct solutions 
at zero temperature that can be associated with the a and 7 phases of cerium. However, despite 
the relatively good structural and magnetic properties, PBEO predicts the 7 phase to be the stable 
phase at ambient pressure and zero temperature, in contradiction with low temperature experiments. 
EX+cRPA reverses the energetic ordering, which emphasizes the importance of correlation for rare- 
earth systems. 



The hrst-principlcs description of / electron materials 
is a considerable challenge and a highly debated topic 
in condensed-matter physics. The simultaneous presence 
of itinerant spd states and localized partially occupied / 
states and their mutual interaction in rare-earth materi- 
als give rise to a rich variety of physical phenomena that 
continue to be a testing ground for electronic-structure 
theories. Cerium is one of the most prominent represen- 
tatives in this regard and, even more intriguingly, under- 
goes an isostructural (fee) a-7 phase transition accompa- 
nied by a volume collapse of 15% at room temperature 
and ambient pressure [J 13. The a phase is character- 
ized by enhanced Pauli paramagnetism and has a smaller 
volume, while the larger-volume 7 phase follows a Curie- 
Weiss behavior for the magnetic susceptibility. 

At zero temperature, first-principles calculations have 
so far been unable to produce a double minimum in 
the total energy versus volume curve, that would be a 
direct indication of the phase transition, within a sin- 
gle theoretical framework. In local or semilocal (LDA 
or GGA) functionals of density-functional theory (DFT) 
the / electrons are always delocalized, due to the strong 
self-interaction error of the functionals, and only the a 
phase is described with some confidence 0, H| . The self- 
interaction corrected local spin density approximation 
0, IE 0] and Hubbard U augmented local or semilocal 
DFT calculations (LDA/GGA+U) 0,0] enforce localiza- 
tion of the / electrons. They subsequently yield a phase, 
whose volume and magnetic moment are consistent with 
the 7 phase, but the description of the a phase requires 
a different treatment, namely LDA. Dynamical mean 
field theory in combination with LDA (LDA+DMFT) 
has been applied to the study of the phase transition at 
finite temperatures fjl-[ll|. but could so far not be ex- 
tended to the zero temperature limit. Whether a double 
minimum exists or would only emerge in the free energy 
curve at finite temperature, due to entropic effects as sug- 
gested by Amadon et al. |lfjj, is therefore still a matter 
of debate. 



In this Letter we show that hybrid density function- 
als 12Hl4j|. that incorporate a fraction of exact ex- 



change yield a double minimum within a single theo- 
retical and computational framework. The results are 
further improved quantitatively by employing exact ex- 
change plus correlation in the random-phase approxima- 
tion (EX+cRPA) (see Ren et al. [PoJ and references 
therein) [3] . In our approach all electrons are treated 
on the same quantum mechanical level, in contrast to 
LDA/GGA+J7 or LDA+DMFT studies. We obtain two 
distinctly different solutions, whose structural, electronic 
and magnetic properties are consistent with experimental 
results for the a and 7 phase, respectively. 

All calculations in this work were performed with the 
all-electron code FHI-AIMS (Fritz- Haber-Institut ab initio 
molecular simulation) [It], [HI , that is based on numeric 
atom-centered orbitals. Relativistic effects are treated at 
the level of the scaled zero-order regular approximation 
fl~9| . Here we present results obtained using the PBEO hy- 
brid functional [13| for both cluster and periodic systems 
[20j and show that the HSE hybrid functional [l4[ yields a 
similar description. For comparison we also applied the 
local-density approximation in the parameterization of 
Perdew and Zungcr 21 1 and the Pcrdcw-Burke-Ernzcrhof 
generalized gradient approximation (PBE) [22j ■ Periodic 
calculations were performed with a 6 x 6 x 6 k mesh. This 
gives energies that are converged to within 5 meV, which 
is sufficient for the energy scale of interest here [cf. Fig. [I] 
(c)]. The hybrid functional calculations were carried out 
with a tier 1 numeric atom-centered orbitals basis [l7| . 
whereas for (EX+cRPA) @PBE0 it proved necessary to 
go up to tier 3 [23| . In both cases ferromagnetic ordering 
is assumed in our spin-polarized calculations. 

In cerium the 4/, 5d, and 6s states all lie in the vicinity 
of the Fermi energy, giving rise to different electronic con- 
figurations that are very close in energy. In cases like this, 
approaches that are based on the density matrix rather 
than the density, such as DFT+J7, hybrid functionals 
or Hartree-Fock are more susceptible to local minima in 
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FIG. 1. Cohesive energy [E coh = —{E - £-E atom )] of cerium 
as a function of the lattice constant (ao). Dashed lines show 
HSE06 results. The spin moment increases with volume for 
the LDA and PBE solutions, while in PBEO and HSE06 it 
remains approximately constant to zero and one-half for the a 
and 7 phases, respectively. Experimental lattice parameters 
for the two phases at finite temperature [I| are marked by 
black arrows. 



the potential-energy landscape of the electrons [24H30j . 
This fact can be exploited to search for distinct, stable 
electronic configurations. For a given lattice constant 
ao our PBEO calculations are initialized with the den- 
sity matrix from a preceding PBE calculation and with a 
high electronic temperature (T ^1000 K). The electronic 
temperature yields a broadening of the one-electron en- 
ergy levels. In subsequent calculations the temperature is 
gradually reduced until a particular solution is stabilized 
at T = 0. When scanning a range of lattice constants 
(ao) for cerium metal, the binding energies from differ- 
ent values of ao fall into two different smooth curves. In 
practical calculations, the solution at one particular lat- 
tice constant can be used to initialize the calculations at 
a neighboring lattice constant. In this way, the two bind- 
ing energy curves can be stabilized with relative ease. 

It is one of the core results of this work that the above 
technique provides two different stable solutions. In Fig. 
Q]the cohesive energies (.Ecoh) obtained from LDA, PBE, 
PBEO, and HSE06 are presented as a function of the 
lattice constant. Our LDA and PBE results are in agree- 
ment with previous calculations 0, H| and exhibit only 
one minimum. The associated volume is consistent with 
the a phase, although the actual value is underestimated. 
Constraining the magnetic moment does not introduce a 
second minimum. In contrast, in PBEO and HSE06 two 
stable solutions are found. One solution has a minimum 
approximately coinciding with the LDA or PBE mini- 
mum, while the second assumes its minimum at a much 
larger lattice constant, consistent with the one of the 7 
phase. The magnitude of the cohesive energy systemati- 
cally reduces from LDA to PBE, and to PBEO. 



The two PBEO solutions differ in their electronic struc- 
ture as e.g. reflected in the density of states (see Rcf. 
[3l| ) and the magnetic moment m (cf Fig. [3] (c)). m of 
the low volume phase lies around 0.2 fj, , while in the high 
volume phase m is close to one. A rapid change of the 
local magnetic moment across the a-7 phase transition 
was also observed in LDA+DMFT [|| . Also the number 
of / electrons is approximately one in both phases, as 
suggested by positron annihilation experiments [32| . Fig- 
ure [2] shows the difference of the PBEO electron density 
between the a- and the 7- like solutions, n a (r) — n 7 (r), 
projected onto a volume slice at the same lattice con- 
stant of 4.6 A. The a-like phase has a higher density in 
the interstitial region, a clear indication of electron der- 
ealization, whereas in the 7-phase solution the electron is 
more localized around the Ce atom. More importantly, 
however, the density difference has the shape of an / or- 
bital of xyz, z{x 2 — y 2 ) symmetry, as evidenced by its 
three-dimensional projection (not shown). This provides 
a strong indication that the balance between localization 
and delocalization of the / electrons plays a key role in 
the emergence of the double minimum in the cohesive 
energy curve. 

Further inspection of Fig. [1] reveals that the energetic 
ordering of the two PBEO solutions is not consistent with 
the experimental phase diagram [33|, which shows that 
the a phase is lower in energy than the 7 phase at low 
temperatures. The opposite is true for the PBEO results. 
To overcome this discrepancy we resort to a more ac- 
curate treatment of exchange and correlation. As the 
name indicates, within EX+cRPA the exchange term is 
treated exactly (i.e. not reduced by a factor as in hy- 
brid functionals) , and correlation is treated at the level 
of the random-phase approximation. The mixing factor 
in the hybrid functionals that controls the fraction of ex- 
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FIG. 2. Volume slice through the difference between the bulk 
Ce electron densities of the a and 7 phases at the same lattice 
constant of 4.6 A, at which both phases have the same energy 
(see Fig. [T]). The a phase has a larger contribution in the 
interstitial region, whereas the 7 phase is more localized and 
exhibits a clear /-orbital shape. 
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FIG. 3. PBE0 results for cerium clusters and corresponding 
bulk values for (a) the cohesive energy (£ co h), (b) the lattice 
constant (ao) and (c) the magnetic moment on the central 
atom (m). All clusters exhibit two solutions that converge 
to the calculated bulk limit. Experimental results marked on 
the right axis are taken from Ref. 



act exchange may be regarded as a simplified screening 
function, which is replaced by a physical and system- 
dependent screening in EX+cRPA. EX+cRPA is so far 
only implemented for finite systems in the current version 
of FHI-AIMS. We therefore adopt the strategy of modeling 
bulk Ce with systematically increasing cerium clusters. 
The clusters are cut from the fee crystal structure, with 
one atom in the center surrounded by shells of first, sec- 
ond, and third nearest neighbors (i.e. thirteen, nineteen, 
and forty-three atoms in total). To reduce edge effects, 
we use the formula for evaluating the effective cohesive 
energy of clusters (see, e.g., [34l lijq) 



E C oh — — 



E ]T(7V c £ c a 




(1) 



where E is the total energy, N c the number of atoms 
in the cluster with c nearest neighbors, and E^ tom is 
the atomic total energy for a c-fold-coordinated atom. 
Basis-set superposition errors are corrected when evalu- 
ating E* tom . For each cluster two distinct solutions were 
always found. Figure [3] demonstrates how the two sets 
of PBE0 cohesive energies, equilibrium lattice constants, 
and magnetic moments (on the central atom) of the clus- 
ters converge towards the corresponding bulk values as 



the cluster size grows. It can be seen from Fig. [3] that 
the essential physics is already captured with the small- 
est cluster size. In the following, the (EX+cRPA) @PBE0 
results for the Cei9 cluster will be shown, which suffices 
for the discussion below. 

The EX+cRPA calculations arc performed non-self- 
consistently with the orbitals and eigenenergies of the 
two PBE0 solutions as input. Figure 2] shows the 
(EX+cRPA) @PBE0 cohesive energy for the 19-atom 
cluster. Since the density of states near the Fermi level is 
higher in the low-volume phase, and consequently screen- 
ing and the RPA correlation energy are larger, the low- 
volume phase moves down in energy relative to the high- 
volume one and the correct energetic order is restored. 
According to the extrapolation of the experimental data 
13, 33 1, the difference in internal energy (At/) between 
the two phases should lie between 20 and 30 meV. 
The (EX+cRPA) @PBE0 value for the 19-atom cluster 
amounts to At/ ~ 45 meV. Although the difference is 
larger than the experimental estimation of the maximum 
energy difference between the two phases, it is compara- 
ble to the entropy contribution TAS at ambient condi- 
tions 33j. Therefore, our results on the electronic con- 
tribution to the phase transition do not rule out that en- 
tropy might play a noticeable role in the phase transition 
at finite temperature as proposed in Ref. Hoi . The com- 
mon tangent construction to the (EX+cRPA) @PBE0 co- 
hesive energy curves leads to a transition pressure of 
Pt — —0.74 GPa at zero temperature, in good agreement 
with the extrapolated experimental Pt ~ —0.8 GPa. The 
calculated lattice constants for the a- and 7- like phases 
are 4.45 and 5.03 A, respectively. This corresponds to a 
volume collapse of ~ 30% at zero temperature, to be con- 
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FIG. 4. Calculated (EX+cRPA) @PBE0 cohesive energy 
(-Ecoh) for the 19-atom fec-cerium cluster as a function of 
ao. The dashed line illustrates the Gibbs construction for 
the transition pressure in good agreement with the extrap- 
olated experimental Pt — —0.8 GPa (HU). Arrows on the 
energy axes: experimental cohesive energy from [3^ ]. 
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trastcd with the 15% observed experimentally at ambient 
conditions [![. 

Over the years many experimental and theoretical 
studies have addressed the origin of the transition but 
a conclusive answer is still lacking. The two prevalent 
propositions are: a Mott transition for the f electrons 
37 . I^H and the Kondo volume collapse [23, Hjj. In 
the Mott picture the hybridization between / orbitals 
is believed to change across the transition, leading to 
one phase with delocalized / electrons (a phase) and the 
other with localized / states (7 phase) . In the Kondo vol- 
ume collapse model, the / electrons are assumed to be 
localized in both phases, and the change in spin screen- 
ing of the localized moments by the conduction spd elec- 
trons is responsible for a change in the Kondo tempera- 
ture across two orders of magnitude, with an associated 
change in system properties. Although the two mod- 
els are different in nature, recent theoretical works sug- 
gest that the resulting scenarios are quite similar at finite 
temperature [4lT | and both are consistent with available 
experimental data [H, EH ■ The common belief that the 
phase transition should be driven by changes in the elec- 
tronic structure alone has also been questioned, but the 
contribution of e.g. lattice vibrations has been estimated 
to be lower than 30% of the total energy change [43M45| . 
Since our calculations are performed at zero temperature 
they can only describe a pressure induced phase transi- 
tion at zero temperature but not a temperature induced 
transition at ambient pressure as observed in the experi- 
mental phase diagram. However, the occurrence of a dou- 
ble minimum in the cohesive energy curve and the strong 
signature of / states in the accompanying density differ- 
ence shown in Fig. [2] suggest that the localization of / 
electrons is a strong contributing factor to the phase tran- 
sition. Density-functional theory (being a ground state 
theory) does not need to give the right electronic excita- 
tion spectrum to describe the structural phase transition 
at zero temperature, even if the spectrum is essential 
for the distinction between the Mott transition and the 
Kondo volume collapse at finite temperature. Neverthe- 
less, our zero temperature results are consistent with the 
scenario of a Mott transition and show that advanced 
DFT exchange-correlation functionals can indeed capture 
such Mott physics. 

In conclusion, PBEO hybrid functional calculations 
combined with EX+cRPA produce multiple distinct so- 
lutions of the electronic structure of bulk cerium at zero 
temperature. These can be discriminated by their mag- 
netic moment and the degree of / electron localization 
and have been tentatively associated with the a and 7 
phases of cerium. At the PBEO level, the energetic order- 
ing of the two solutions is reversed compared to the zero 
temperature extrapolation of the experimental phase di- 
agram. EX+cRPA recovers the right ordering, which 
highlights the role of correlation in rare-earth systems. 

An interesting aspect that emerged from the cluster 



extrapolation approach is the presence of a volume col- 
lapse at the nanoscale down to the dimer. This is a 
very unusual feature for a first-order phase transition in 
a metal and opens the question of whether other lan- 
thanides would exhibit the same behavior. This predic- 
tion suggests further investigation both theoretically and 
exp er iment ally. 
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